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New classes of topological crystalline insulators with unpinned surface Dirac cones 
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We theoretically predict two new classes of three-dimensional topological crystalline insulators 
(TCIs), which have an odd number of unpinned surface Dirac cones protected by crystal symmetries. 

The first class is protected by a single glide plane symmetry; the second class is protected by a 
composition of a twofold rotation and time-reversal symmetry. Both classes of TCIs are characterized 
by a quantized tt Berry phase associated with surface states and a Z 2 topological invariant associated 
with the bulk bands. In the presence of disorder, these TCI surface states are protected against 
localization by the average crystal symmetries, and exhibit critical conductivity in the universality 
class of the quantum Hall plateau transition. These new TCIs exist in time-reversal-breaking systems 
with or without spin-orbital coupling, and their material realizations are discussed. 


The notion of symmetry protected topological (SPT) 
phases has recently emerged from studies on topolo, 


og- 

S- 


ical insulators and is now being intensively studied^ 

. SPT phases generically have gapless boundary states 
that are stable against perturbations, provided that 
certain symmetry is preserved. The topological prop¬ 
erty of these boundary states depends crucially on the 
underlying symmetry. In the well-known example of 
three-dimensional (3D) topological insulators, time re¬ 
versal symmetry protects an odd number of surface Dirac 
points that are pinned to time-reversal invariant mo¬ 
menta (TRIM). 

Recent theoretical studies [H-flsj] have found a variety 
of SPT phases that are protected by crystal symme¬ 
tries, termed topological crystalline insulators (TCIs)0. 
A universal property of TCI phases is the presence of 
protected surface states on symmetry-preserving sur¬ 
faces. However, depending on the underlying crystal 
symmetry, surface state properties of different classes 
of TCIs can vary significantly. In spin-rotationally- 
invariant TCIs protected by rotational symmetries of 
the crystal (such as Q)!, [H , surface states exhibit 
band crossings that are pinned to certain high-symmetry 
point(s) in the SBZ. In spin-orbit-coupled TCIs protected 
by mirror symmetry as realized in the SnTe class of IV- 
VI semiconductors!^ I3“22|, surface states exhibit Dirac 
points on a specific mirror-symmetric line in the SBZ, 
corresponding to the projection of the 2D plane with a 
nonzero mirror Chern number in the Brillouin zonej^. 

In this work, we theoretically predict two new classes 
of 3D TCIs, protected by a glide plane symmetry and 
a space-time inversion symmetry respectively. Unlike all 
topological insulating phases known so far, their surface 
states consist of a single (more generally an odd number), 
unpinned Dirac point with a quantized 7r-Berry phase. 
Importantly, these new TCI phases are robust against ei¬ 
ther magnetic or nonmagnetic impurities, which by def¬ 
inition preserve the crystal symmetry on average. Re¬ 
markably, the disordered surface realizes, without any 
tuning, a critical phase in the universality class of quan¬ 
tum Hall plateau transition. 

TCI with glide symmetry The first class of TCI exists 


in 3D systems (with or without spin-orbital coupling) 
that have a glide plane symmetry, i.e., a combination of 
reflection and a translation by half a lattice vector: 

Mg'. {x,y,z) ^ [x,y,-z)+& 1 / 2 , (1) 

where {x, y, z) is the position vector and B. 1 ^ 2,3 are the 
basis of lattice vectors, out of which ai is inside the xy- 
plane. A key difference between a mirror plane and a 
glide plane is that the mirror plane squares to identity 
(up to a Berry phase associated with a 27r rotation), while 
the glide mirror squares to a lattice translation: 

M2=(-1)/Ta,. (2) 

Here / = 0 applies to spin-rotationally-invariant systems, 
where reflection does not involve spin; / = 1 applies to 
spin-orbit-coupled systems, where reflection acts on spin 
s = i and squares to —1. 

We now show how a glide plane can protect a crossing 
point in the surface bands. First, a symmetry-preserving 
surface must be (i) perpendicular to the glide plane {xy- 
plane in this case) and (ii) invariant under the translation 
along ai. Without loss of generality, we choose ai ^2 to 
be along x, y-axes, respectively. The only surface that 
satisfies both conditions is then the x^-plane, whose SBZ 
is plotted in Fig. [Ija). Due to the translational sym¬ 
metry in the xz-plane, the Hamiltonian with an open 
surface is diagonal in Bloch basis with crystal momen¬ 
tum (kxykz). The corresponding Bloch Hamiltonian is 
denoted by h{kx,kz). The presence of glide plane sym¬ 
metry implies: 

MG{kx)h{kx,kz)MG^ {kx) = h{kx,-kz). (3) 

Here the operator Maikx) represents the action of Mq 
in Bloch basis. Note that unlike point group symme¬ 
try operators, AIg(^x) is a function of kx- This results 
from Eq.(l2]), which implies eigenvalues of the glide plane 
are not constants but depend continuously on electron’s 
momentum. 

It follows from Eq.(l3|) that all bands on the two high- 
symmetry lines fcz = 0 and tt (where lattice constants are 
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taken to be unity) can be labeled by the eigenvalues of 
■Mc{kx)- Using Eq.Q, and taking into account that at 
7ai = the two eigenvalues of Mc{kx) are 

m±{kx) = ( 4 ) 

These eigenvalues divide the bands along k^ = 0,7r into 
two branches that have glide plane eigenvalue of m+(kx) 
and m-[kx) respectively—hereafter referred to as m_|_ 
and m_ bands. Since any hybridization between a m_|_ 
band and a m_ band breaks the glide plane, a single 
crossing point ko between them at any momentum on 
the line = 0 or = tt is protected. The k ■ p Hamilto¬ 
nian at ko takes the form (up to a unitary basis change): 
H(k) = VxkxUi + VykyCr 2 -, where k denotes the momen¬ 
tum relative to the Dirac point; and the action of glide 
mirror on the two degenerate surface states at kp is rep¬ 
resented by McCko) cx CTi, up to a C/(l) phase factor. 
Provided that glide plane symmetry is preserved, pertur¬ 
bations can shift the band crossing point kp along the 
high-symmetry line, but cannot open a gap. This leads 
to a symmetry-protected surface Dirac point that is not 
pinned to a specific point. Furthermore, the stability of 
this Dirac point has a topological origin arising from the 
quantization of the Berry’s phase, a point to which we 
will return later. 

Having addressed the protection of a single surface 
Dirac cone, we proceed to study the stability of two Dirac 
cones, each of which is locally protected by Mq- In this 
case, can we adiabatically tune the Hamiltonian to fully 
gap the surface without closing the bulk gap? This turns 
out to be a subtle issue that has not been encountered in 
topological insulating phases studied so far. To answer 
this question, two cases should be discussed separately: 
(i) the two Dirac points appear on the same line, either 
kz = 0 or kz = TT and (ii) there is one Dirac point on each 
mirror invariant line. 

As an example of the first case, we combine two identi¬ 
cal copies of the aforementioned TCI with a single surface 
Dirac cone. The combined system then has two surface 
band crossings that appear at the same point on same 
mirror symmetric line in SBZ, say (fco,0) [Fig. [Ijb)]. In¬ 
finitesimal perturbation on the surface cannot fully gap 
the spectrum, because the two right-going (left-going) 
modes have the same mirror eigenvalue. Two Dirac 
points are hence ‘locally stable’. However, we find the 
corresponding surface state spectrum can be and can only 
be gapped out by sufficiently strong deformations, i.e., it 
is globally unstable. To see this, we need to use a key 
property of the glide plane symmetry, as indicated in 
Eq.Q: a m_|_-band is connected with a m_-band at the 
BZ boundary at kx = ±7r, because the phase factor on 
the right hand side gives an additional minus sign when 
kx goes to fca; -|-27r. Consider a finite surface perturbation 
that pushes ki to the left and k2 right [FigHDc)]. When 
they meet each other again at the SBZ boundary, ac¬ 
cording to the above property, the right-going (left-going) 


modes have opposite mirror eigenvalues so fci _2 can anni¬ 
hilate each other [Figjljd)]. We emphasize that both the 
local stability and global instability of two Dirac points 
are key characteristics associated with the glide plane 
symmetry, which have not appeared elsewhere. We fur¬ 
ther note that since two Dirac points can only annihilate 
each other by crossing the SBZ boundary, the new TCI 
phase with glide mirror symmetry cannot be treated in 
continuum models where fc-space is effectively a sphere 
rather than a torus. Generally, the torus nature of SBZ 
must be considered in studying nonsymmorphic symme¬ 
tries of a lattice. 

In the second case, we consider the spectral flow 
of the band dispersion in the SBZ along the path 
X'TXMZM'X', shown in Fig. (Ha) by the arrows. 
Fig. |2Ka) and FiglJKc) show typical spectral flows for 
a trivial and a nontrivial phase, respectively. We need 
the following principle for the analysis: along X'TX and 
M' ZM, the bands must appear in pairs that cross each 
other. The proof of the principle is give n in Sec. I of 
Supplementary Materials (also see Ref.[23j). In Fig.|2Ka), 
there are two band crossings on X'TX and M'TM, re¬ 
spectively. When the chemical potential on the surface 
increases, the two bands move upward in energy together, 
as in Fig. |2](b) and are eventually pushed into the con¬ 
duction bands as in Fig. He), leaving a full gap on the 
surface. In Fig. HKd), there is one band crossing along 
X'TX, and no crossing along M'ZM. When the surface 
chemical potential increases, the two bands move upward 
in energy together. However, since bands must appear in 
crossing pairs along kz = 0 and kz = tt, pushing up the 
chemical potential will ‘pull out’ a pair of bands from the 
valence bands, as shown in Fig. He). As as result, the 
Dirac point on X'TX moves to M'ZM without closing 
the bulk gap, and the flow remains [see Fig.[2][f)]. 

We have now shown that a single surface Dirac cone 
can be protected by a glide mirror symmetry, but not 
two cones. Due to the bulk-edge correspondence, this 
implies the existence of a Z 2 topological invariant for 
the bulk. The analytic and explicit expression of the in¬ 
variant involve the theoretical tool of non-Abelian Berry 
phase and Wilson loop which we leave to Sec. B 

of Supplementary Materials. There, we also provide a 
lattice model as an example of this TCI phase in a spin¬ 
less system. 

To summarize, the TCI with a glide mirror symmetry 
has a single Dirac point that (i) is unpinned to any high- 
symmetry point, and (ii) under symmetry perturbations 
that does not close the bulk gap, can move along two 
mirror invariant lines as well as shift between the two 
lines. Can we have a new TCI phase with a single Dirac 
point that can freely move in the SBZ without being 
confined to any high-symmetry point or line in the SBZ? 
Below we provide an affirmative answer to this question. 

TCI with space-time inversion symmetry We now show 
that an anti-unitary symmetry C 2 *T, the combination of 
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FIG. 1: (a) The SBZ of the xz-surface of a 3D (spinless) system with glide plane defined in Eq.([T]). (b,c,d) are schematics 
of the process where two identical Dirac points on a single mirror symmetric line annihilate each other through finite surface 
perturbation. Here in (b) we set fci ,2 slightly away from each other only to indicate that there are two instead of one Dirac 
points. The color bar on the left of (b) shows the color code for the phase of the glide plane eigenvalues. 
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FIG. 2: (a-c) show the process of gapping the spectral flow 
when there are two band crossings, on X'TX and M'VM 
respectively, (d-f) show the continuous flow when there is 
only one band crossing, and how the crossing moves from 
X'fX to M'VM. 


a twofold rotation and time-reversal, protects a Z 2 TCI 
phase with a single Dirac cone on the surface perpen¬ 
dicular to the twofold axis, whose location in the SBZ 
is completely unpinned. This symmetry is hereafter re¬ 
ferred to as a space-time inversion symmetry, ds C 2 *T 
operation on the surface sends (x, y, t) to (—x, —y, —t). 

We first note the identity (C 2 * T)^ = /, which holds 
for both systems with and without spin-orbit coupling, 
because both Cf and are equal to —I for the former, 
and / for the latter. This identity allows us to repre¬ 
sent the action of C 2 *T on the two degenerate states 
at a (putative) Dirac point by C 2 * T = Ffcri, where K 
is complex conjugation. Then, one can easily verify that 
the Dirac Hamiltonian H(k) = Vxkx(^i+Vykya2 is indeed 
compatible with C 2 *T symmetry. Perturbations that are 
invariant under C 2 *T correspond to cti or a 2 , which sim¬ 
ply shift the Dirac point. The Dirac mass term (T 3 is odd 
under C 2 *T and hence forbidden. This simple analysis 
demonstrates the protection of an unpinned Dirac point 
by C 2 * T. 

We now further show that two cones can annihilate 
each other under perturbations preserving the C 2 * T 
symmetry. The Hamiltonian for two identical copies of 
surface Dirac cones has the following form H = tq 0 


{kx<Jx + kyay), with C 2 *T = ifroGicri, where 70 , 1 , 2,3 are 
identity and Pauli matrices acting on the species space. 
A symmetry preserving perturbation 5 mT2 1 X 3 can gap 
the whole spectrum. This implies the Z 2 nature of the 
TCIs protected by (72 * T symmetry, with an odd number 
of surface Dirac points. The Z 2 topological classification 
can also be derived from analysis of the bulk topology, 
as we show in Sec. HI of Supplementary Materials. 

The protected surface states of TCIs with space-time 
inversion symmetry are characterized by a quantized 
TT Berry phase, for any loop enclosing all Dirac points 
thereof. As stated above, the same is true for TCIs with 
the glide mirror symmetry, provided that the loop en¬ 
closing the Dirac point is symmetric. Here we provide 
a general proof for both cases. First we show that the 
Berry’s phase is independent of the geometry (such as 
size or shape) of the loop. According to Stokes the¬ 
orem, the difference between Berry phases associated 
with two loops mod 27r is given by the integral of Berry 
curvature over the region bounded by them. For sys¬ 
tems with glide plane symmetry, the Berry’s curvature 
for the surface states satisfies F{kx,kz) = —F{kx,—kz), 
due to the oddness of F under reflection z —>• —z. So 
for any loop symmetric about a mirror invariant line, 
the Berry’s curvature integral inside the loop subtract 
the singular band crossing point vanishes. For systems 
with space-time inversion symmetry. Berry’s curvature 
satisfies F{kx,ky) = —F{kx,ky) due to the oddness of 
F under time-reversal symmetry, and hence vanishes ev¬ 
erywhere in k space. So the Berry’s curvature integral is 
again zero. Therefore, only Dirac points make a singular 
contribution to the Berry phase. As is well-known, the 
Berry phase of each Dirac point is tt. Therefore, an odd 
number of Dirac points on the TCI surface guarantees a 
quantized tt Berry’s phase, for both glide mirror or space- 
time inversion symmetry. (The quantization of tt Berry 
phase hy C2 *T and its implication for TCI phase were 
mentioned in a recent work [ 2 ^; The local stability of a 
Dirac point in the presence of C2 *T symmetry was also 
noted 3^, 311.) 

Material realizations The two new classes of Z 2 TCI 
can be realized in systems with or without spin-orbit cou- 
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pling. They are discussed separately below. 

In spinful systems, the two new classes of TCI are 
consistent with time-reversal symmetry (TRS). One can 
simply consider a Z 2 strong topological insulator which 
also has glide mirror and/or twofold rotation symme¬ 
try. Pick a surface that preserves the glide plane or the 
twofold axis. It is guaranteed to possess a single Dirac 
cone at one TRIM, protected by spinful TRS. Now let 
us add perturbations that break TRS but preserve the 
glide plane or C 2 * T; we immediately obtain the Z-i TCI 
phases found in this work. This observation makes it very 
simple to find these new TCIs in spinful systems. For 
example, the newly discovered topological Kondo insula¬ 
tor SmBfi fs^ - l^ has C' 2 -symmetry on the (OOl)-surface. 
Adding any magnetic field parallel to the surface pre¬ 
serves C 2 * T. Hence SmBe with an in-plane field can be 
considered as a Z 2 TCI protected by C 2 * T, having an 
odd number of Dirac cones located away from any TRIM. 

To realize the new TCI phases in systems without spin- 
orbit coupling (or equivalently, spinless systems) requires 
breaking TRS. This can be shown by proving that a sin¬ 
gle Dirac cone characteristic of these TCIs cannot ap¬ 
pear on the surface of a time-reversal-invariant spinless 
insulator. To see this, first note that Dirac cones at 
non-TRIM must appear in pairs with opposite momenta. 
This leaves the possibility of having a single Dirac cone 
at one TRIM. However, given that time-reversal symme¬ 
try is represented hy T = K (up to a gauge) for spinless 
fermions and it reverses k measured from the TRIM, a 
2D Dirac Hamiltonian is simply not allowed by T. This 
is because only one of the three Pauli matrices, namely, 
cr 2 , is reversed under K, while a 2D Dirac Hamiltonian 
such as H(k) = VxkxO'i + Vykya 2 or any other form must 
involve two Pauli matrices. This concludes that a sin¬ 
gle Dirac cone cannot exist on the surface of a spinless 
system with time-reversal symmetry. 

Therefore, the key requirement for spinless TCIs is to 
break TRS while preserving the relevant crystal symme¬ 
try of either glide mirror or space-time inversion. One 
may search in magnetic insulators with negligible SOC, 
since many types of magnetic order have at least one of 
the two symmetries. Alternatively, we note that photonic 
crystals may be a very promising platform for finding 
these TCI phases, because their structures and crystal 
symmetries can be easily manipulated [s^, 

Stability against disorder and quantum Hall criticality 
We now show that topological surface states of TCIs with 
either glide mirror or (72 * T symmetry are fully robust 
against any type of disorder (magnetic or nonmagnetic) 
and cannot be exponentially localized even under strong 
disorder on the surface. Similar to the case of a disor¬ 
dered TI surface with random magnetic impurities(38|. 
the delocalization of TCI surface states here is protected 
by the average symmetry!^ [H, ” by definition the 

ensemble of TCIs with all disorder realizations must pre¬ 
serve the relevant symmetry. Consider a finite area on the 


surface where surface states are gapped by disorder. The 

mass m acquired by the Dirac cone can either be positive 

or negative, corresponding to two types of domains. The 

symmetry (glide mirror or C 2 *T) maps am > 0 domain 

to a m < 0 domain and vice versa. Importantly, the two 

2 

domains have quantized half-integer Hall conductivity ^ 
2 

and respectively; between the two domains of oppo¬ 
site signs there exists a chiral edge mode. Therefore, in 
the presence of random disorder, these chiral modes per¬ 
colate through the surface and the surface states remain 
delocalized. This argument also shows that the disor¬ 
dered TCI surface is symmetry-enforced to be in a critical 
phase (rather than fine-tuned to a critical point) in the 
same universality class as the quantum Hall plateau tran¬ 
sition. Therefore, the surface exhibits a universal longi- 

2 

tudinal conductivity on the order of ^, and a nontrivial 
scaling of longitudinal and Hall conductivity as a function 
of an applied out-of-plane magnetic field B, which gaps 
the Dirac point and drives each surface into a quantum 
Hall state with a Hall conductance sign(i?)e^/2/i. The 
space-time inversion symmetry protects the gapless sur¬ 
face Dirac point by forbidding this out-of-plane magnetic 
field, thereby defining a TCI phase. Other symmetries 
such as reflection or any improper rotation can also for¬ 
bid the perpendicular field, and therefore guarantees at 
least a Z 2 classification of TCIs0 , as in the two classes 
studied here. 


To conclude, we have theoretically predicted two new 
classes of 3D Z 2 TCI that have unpinned surface Dirac 
cones, which are protected by a glide plane and C 2 * T, 
respectively. The Z 2 nature distinguishes the new TCI 
from the Z TCI protected by mirror symmetry such as 
SnTe. Because of the ‘unpinned’ nature of the surface 
Dirac cone, the TCI phases studied in this work lie be¬ 
yond the previous classification of TCIs with nonsymmor- 
phic space groups based solely on considerations of 2D 
irreducible representations of the wallpaper groups [14 1. 
Besides looking at the surface states, we also mathemat¬ 
ically prove the Z 2 classification which directly reveals 
the bulk band topology. We emphasize that the new TCI 
phases can be realized in both spinful and in spinless sys¬ 
tems. In spinful systems, these TCIs can be realized by 
applying perturbations that break TRS yet preserving 
glide plane or C 2 * T. The realization of TCIs in spinless 
systems is an interesting subject that we leave to future 
work. 
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PROOF OF A PRINCIPLE FOR THE BAND CROSSINGS IN THE PRESENCE OF A GLIDE PLANE 

SYMMETRY 

In the text, we mention that along a mirror invariant line in the SBZ, all edge modes must appear in pairs that 
cross each other, in the absence of chiral modes. We use k S [0,7r) to parameterize this line. All edge modes along 
this line can be separated into m+-bands and m_-bands, and in the main text, we have proved that 

E+,{k = 2^) = E_p{k = Q), ( 5 ) 

where E±i{k) denotes the dispersion of the i-th band in the m± sector. If E^iik) does not cross any other bands, 
then due to 


E+,{k = Q)^E+,{k = 2TT), ( 6 ) 

E^i(k) is a chiral mode, contradicting the assumption of the absence of chiral modes. Therefore, there must be one 
other band that crosses E^i(k). But since only a m_-band can cross a TO+-band, E^iik) must cross some E-j{k). 
The principle is proved. 
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FIG. 3: The 3D Brillouin zone and the surface Brillouin zone on the (OlO)-surface. The blue planes are where the mirror Berry 
curvatures are evaluated in Ea. (l7l8ll and the red (half) plane is where the total Berry curvature is evaluated in Eq.®. The 
two Wilson loops are marked with arrows. 


EXPLICIT EXPRESSION FOR THE Z 2 INVARIANT FOR 3D INSULATORS WITH A GLIDE PLAN 

SYMMETRY 

In the main text, we mention that we have fonnd the explicit expression for the Zi invariant using Wilson loops. 
We define two Wilson loops: (I) a straight line from (tt, 0,0) to (tt, 27r,0) and (II) a straight line from (tt, 0,7r) to 
(tt, 27r,7r). (See the two lines denoted by ‘I’ and ‘IF in Figo) On these lines one can separate the occupied spaces 
into two, one spanned by states with mirror eigenvalue +i^~^ and the other with mirror eigenvalue The 

eigenvalues of the Wilson loop, which are gauge independent, are denoted by ■ Then we consider the plane 
defined by = 0 and kz = tt. On these planes, the bands can be diagonalized into blocks having mirror eigenvalue 
m±{kx) = We can define the Berry phases associated with the two blocks to be for the two planes. 

For now we would assume that the Chern number on both planes are zero, so Next we consider a 

plane that has the two loops as its boundary (the red half plane in Fig.Q. The Berry curvature integral on this plane 
is denoted by Bq. 

We can now define three integers. First we observe that on the plane kz = 0, any state at k ^ = —tt continuously 
evolve to a m_ state at k ^ = tt. This is equivalent to saying that all the eigenvalues corresponding to m = +i of loop 
I must evolve into the eigenvalues corresponding to m = —i by going through the kz = 0 plane. At the same time, 
we know that these eigenvalues gain a total phase of B^ in the process. Therefore, we have 

log(A}+)= Y log(A}_)+ 2n77r. (7) 

J = l....,iVo„e/2 i = l,...,7V„ec/2 

Similarly on kz = tt, we can define n//: 

Y ^og{X^Y+Bl^= Y log(A”)+2n//7r. (8) 

j=l,...,N„^cl2 j=l,...,Nocc/'2 

Finally, all eigenvalues of loop I must evolve into all the eigenvalues of loop II by going through the red half plane 
defined above. At the same time, we know that these eigenvalues gain a total phase of Bq. Therefore 

Y [log(A}+) + log(Aj_)] + Bo= Y [log(A}:5_)+log(Aji)]+2no7r. (9) 

j = l,...,Noac/2 j = l,...,Nocc/2 

Finally, we can write down the topological invariant 

Z 2 = (n/ + nij + no) mod 2. (10) 

In order to call this number as a new quantum Z 2 number, we still need to establish that (i) it is gauge invariant, (ii) 
it does not change under gapped, continuous deformation of the Hamiltonian and (iii) there is at least one example 
that realizes the nontrivial phase, i.e., Z 2 = 1. 

(i) It is trivial to see the gauge invariance because it is defined in terms of Berry curvature integrals and Wilson 
loop eigenvalues, both of which are gauge invariant. 
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FIG. 4: The surface band structure for the xz-surface of the system given in Eq. m for the parameter set m = 2.5, to = t'o = 
t'o =t = t' = t" — 1 and (f> = 0.4. 

(ii) Upon a gapped, continuous deformation of the Hamiltonian, the Wilson loop eigenvalues change continuously. 
However, due to the logarithms used in Ea. (j7l8l9|) . we should still prove that the Z 2 is unchanged even if any of the 
eigenvalues crosses the branch cut. We assume that the branch cut is made such that the logarithm takes value with 
[0, 27r). Suppose one of the eigenvalues Aj_|_ increases its phase that it crosses the cut from below, then it is easy to 
see that n/ ^ n/ + 1 and no ^ no + 1 simultaneously, so Z 2 is unchanged. 

(iii) The tight-binding Hamiltonian is given by 

H{k) = {m — tocoskx — t'ocosky — t"oCOskz)T,Q3 + tsm[{kx — (t>)/‘2,][smkx/2I]ii + coskx/ 2 I] 2 i] (H) 

-I- t' sin fcy So 2 -b t" sin S 31 , 

where S^ = Ci ® aj. The glide mirror operator is given by 

Moikx) = e-^’^-/^[cos{kx/2)Ew + sin(A:,/2)S2o]. (12) 

It can be easily verihed that 

M% = ( 13 ) 

Mcikx + 27 t) = M{kx), 

.A^qH (^kx ^ ky ^ kz^.AA Q — ll{kx:ky^ k^)- 

Fig. |4] shows the surface states on the (OlO)-surface for a certain set of parameters (see caption). There is a single 
Dirac cone in the SBZ along X'TX, which is consistent with the nontrivial spectral flow discussed above. This model, 
with a certain set of parameters, realizes a spinless Z 2 TCI protected by glide plane symmetry. 


PROOF OF THE Z2 CLASSIFICATION FOR 3 D INSULATORS WITH C2 * T FROM BULK TOPOLOGY 

In any tight-binding model, the matrix representation of C 2 * T is KU, where U is some unitary matrix. Since 
(C 2 *T)'^ = I for both spinless and spinful fermions, U must also be symmetric, thus being orthogonal. We can hence 
choose another basis, which is related to the original basis by matrix O, such that O^UO = I. In this new basis, 
C 2 is complex conjugation, or, C 2 *T = K. C 2 *T then ensures that H{kx, ky,kz) = H*{kx, ky,—kz), where H{k) 
is a 3D gapped Hamiltonian in the same basis. Here the sign before kz is inverted hy C 2 * T because C 2 does not 
change kz and T flips its sign. On kz = 0 and kz = planes, the Hamiltonian is real. At each point on these planes, 
the real Hamiltonian along with the Fermi energy defines a projection from all n + m bands (with real wavefunctions) 
to n occupied bands, or equivalently, choosing n independent vectors from an n -I- m-dimensional real vector space. 
Mathematically, the space of these projections is called real Grassmanian manifold, G^{m,m -I- n), the topology of 
which is well-known in algebraic topology [d^. Specifically, its second homotopy group 7 r 2 [G^(TO,TO-|-n)] = Z 2 , which 
means the = 0, tt planes can be each associated with a Z 2 topological invariant. These two Z 2 indices give a Z 2 X Z 2 








classification of the 3D bulk. How did we miss one Z 2 index in the analysis of the surface states? This is because on 
the surface is not a good quantum number, so the two indices cannot be probed distinctly on the surface. We claim 
that the strong index, which indicates the presence/absence of the surface Dirac point is given by the sum of the 
two Z2 indices from kz = 0 plane and k^ = tt plane = S(kz = 0 ) + S{kz = tt), where 6 {kz = O/tt) is the Z2 index. 
The argument goes as follows: if the two planes have the same Z 2 indices, the 3D insulator is adiabatically equivalent 
to layers of decoupled 2D insulators, and cannot have gapless surface modes on the top layer which is insulating. 

We can write down a simple spinless tight-binding model that shows a nontrivial Z 2 phase protected hy C 2 * T: 

H{k) = [m — cos{kx — a) — cos{ky — b) — cos/c^JEga 

-I- sm{kx — n)Sii -f sin{ky — 6 )E 3 i -|- sin (14) 

where 2 < to < 3 and a, b are arbitrary numbers. It is straightforward to check that the Hamiltonian is invariant 
under C 2 *T = K, i.e.. 


H{kx,ky,kz) = H*{kx,ky,kz). 


( 15 ) 


The model of Eg. (11^1) shows a single Dirac cone at (a, b) in the SBZ of the (OOl)-surface. 


